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Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass 

estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. 

Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass 

discrepancy. 

Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core 
overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, 

Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% 
of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in 
the stellar models. 

Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the 
Cepheid mass discrepancy. 
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1. Introduction 

The pulsation properties of Cepheids are tightly correlated to 
their fundamental parameters, such as mass and luminosity, 
which makes them valuable tools for distance measurements 
and cosmology. Cepheids are also powerful probes of stellar 
evolution thanks to the coupling of stellar evolution and stel- 
lar pulsation m odels, both constraining the internal structure of 
these stars (e.g. lHofmeister et dJl"l964UCox ei1rill966llChristvl 
fl966h . However, m ass predictions using each method do not 
agree. IStobiel (Tl969) found that stellar evolution models predict 
Cepheids have higher masses than do stellar pulsation models 
for the same effective temperature and luminosity. This Cepheid 
mass discrepancy has been a challenge for stellar evolution and 
pulsation theory for the past 40 years. 

ICoxl ([1980) showed that the mass discrepancy, defined as the 
mass difference relative to the predicted stella r evolution mass, 
is appro ximately 40% . iMoskalik et all (1 19921) claimed that the 
updated llglesias et al.l (1990) opacities provide a resolution to 
the mass discrepancy. Howev er, the current status of the Cepheid 
mass discrepancy is 17+5% dKeller & Woodl2002l : ICaputo et alJ 
l2003lKeller&W ood 200l lKellerll2008ir Furthermore, there is 
evidence that the m ass discrepancy is a function of both mass 
dCaputo et alJl2005l) and metallicity dKeller & Woodll2006l) . 

Dynamic masses have been determined f or four Galactic 



Cepheids that are in binary systems: SU Cyg (Evans & Bolton 



1990) 



V350 Sgr (lEvans et alJ 
and Polaris (lEvans etaLl 



S Mus (Evans et al 



2006), and Polaris (lEvans et al.l 120081) . The dynamic masses 
are all lower than masses predicted using stellar evolution the- 
ory and consisten t with stellar pulsation models. Furthermore, 
iPietrzvhski et al.l (120101) determine the mass of the Large 
Magellanic Cloud Cepheid OGLE-LMC-CEP0227, which is 
in an eclipsing binary system, to a precision of 1% and find 
that it agrees w i th the mass predicted by stellar pulsation. 
ICassisi & Salarisl d201 ll) model the evolution of this Cepheid and 
find agreement with the dynamic mass when extra mixing is in- 



cluded. These results suggest that there are physics missing in 
the stellar evolution calculations. 

The two most likely solutions to the mass discrepancy are 
convecti ve core overs hooting in a Cepheid's main-sequence pro- 
genitor dChiosi etal.|ll992l) and mass loss during the Cepheid 
stage of evolution (Bon o et al. 2006). Convective core over- 
shooting during main sequence evolution mixes extra hydrogen 
into the core. This leads to a more massive post-main sequence 
helium core, hence to a more luminous Cepheid or conversely to 
a less massive Cepheid for the same luminosity if overshooting is 
not included in the stellar evolution models. Overshooting is also 
required to explain observations of ec lipsing binary stars (e.g. 
ISandberg Lacv et al Il2010t IClausen et al]|2010l) , Cephei stars 
(lLove kin & Goupil 2010) and massive B-type stars (Brott et al. 
120111) . 

On the other hand, mass loss during the Cepheid stage 
of evolution acts to reduce t he ste llar mass without affecting 
the stellar luminosity. iDeasvl (1 19881) determined mass-loss rates 
of 10 ~ 9 to 10~ 8 Mr yr' 1 from IRAS observations. More re- 
cently M erand et al.l (120071 and references therein) observed in- 
frared excess in nearby Galactic Cepheids from interferomet- 
ric observations, while Spitzer observat i ons also detected in - 
frared excesses (|Marengo et al.l l2010albt |Barmbv et al.l 120 lTT) . 
iNeilson et alJ d2009ll2010l) modeled the infrared excess in Large 
Mage llanic Cloud Cepheids in the OG LE-IH dSoszvnski et alJ 

2008) and SAGE dMeixner et al.ll2006l) surveys. In these works, 
the observed infrared excess was modeled by a dusty wind, sug- 
gesting that Cepheids may be undergoing significant mass loss. 

From a theoretical perspective, INeilson & Lesterl (120081 

2009) developed a prescription for pulsation-driven mass loss 
in Cepheids. They predicted mass-loss rates up to 10~ 7 M Q yr~' . 
While this evidence suggests Cepheid mass loss is important, 
it remained unclear whether enough mass is lost during the 
Cepheid stage of stellar evolution to account for the measured 
Cepheid mass discrepancy. 



1 



Neilson et al.: The Cepheid mass discrepancy and pulsation-driven mass loss 



The purpose of this work is to test whether pulsation-driven 
mass loss in Cepheids is an important contributor towards solv- 
ing the Cepheid mass discrepancy. In the next section, we es- 
timate the order-of-magnitude change in mass that may occur 
du ring the Cepheid stage of evolution due to mass loss based on 
the lNeilson & Lesten(l2008l) prescription. In Sect. [3] we compute 
detailed stellar evolution models to explore the role of mass loss 
and convective core overshooting. In Sect. [4] we summarize our 
results. 

2. Analytic test of pulsation-driven mass loss 

We estimate the amount of mass loss during the Cepheid stage of 
evolution, where we assume the Cepheid lifetime is equivalent 
to the helium-burning timescale Te e for a given stellar mass. 

We compute an average Cepheid mass-loss rate by assum- 
ing that the Cepheid instability strip is infinitesimally thin. Thus, 
for a given mass and luminosity there is only one value for the 
pulsation period, amplitude of luminosity variation, and ampli- 
tude of radius variation. The mass-loss rate is then computed 
using the pulsation-driv en mass-loss prescription developed by 
iNeilson & Lester! (12008), who hypothesized that pulsation in the 
envelope of a Cepheid generates shocks that carry momentum to 
the surface of the star and enhances the mass loss that an evolved 
star undergoes. Mass-loss rates were found to be enhanced by up 
to three orders of magnitude. In that prescription, the mass -loss 
rate is a function of the stellar mass, luminosity, radius, pulsa- 
tion period, and pulsation amplitudes of the luminosity and ra- 
dius. We d etermine the pulsation period using the bolometric 
Leavitt law dTurnerl2 010). an d the radius is given by th e period- 
radius relation (iGieren etal.ll989t INeilson et al.l2 010). The am- 
plitudes of the velocity and brigh tness variation are compute d 
from period-amplitude relations (Klagvivik & Szabados 2009). 
The change in radius is the period times the velocity amplitude, 
and we assume the V-band amplitude is equivalent to the bolo- 
metric change of brightness. We note that the predicted pulsation 
amplitudes have significant errors, but the only other way to pre- 
dict these amplitudes is using nonlinear pulsation models. Thus, 
given a stellar mass and luminosity, we determine the period, ra- 
dius, and pulsation amplitudes, hence the mean mass-loss rates. 

The Cepheid's luminosit y is determined fr om stellar evolu- 
tion calculations. We use the Heg er et al.l(l2000l) stellar evolution 
code to compute models with masses M — 4 - 9 M in steps 
of 1 M Q . The models are computed by assuming no convec- 
tive core overshooting in the main sequence progenitors and no 
pulsation-driven mass loss during the post main-sequence evolu- 
tion. We also determine the helium-burning timescales for each 
mass from the models. We use the predicted timescales re e and 
pulsation-driven mass-loss rates M to determine the contribution 
of mass loss towards the mass discrepancy, as shown in Table. Q] 
The range of mass-loss rates is consistent with the results of 
INeilson & Lesterl d2008i) . It should be noted that the lower mass- 
loss rate for the 6 M model is due to a local minimum in the 
period-amplitude relation at a period < 10.4 day. These results 
suggest, however, that mass loss is not a solution for the Cepheid 
mass discrepancy of 17 + 5%. The change in mass relative to 
the initial mass for the low-mass Cepheids appears consistent 
with the discrepancy but for the high-mass (8, 9 M ) Cepheids, 
there is a negligible mass change due to pulsation-driven mass 
loss. Furthermore, the helium-burning timescale is longer than 
the Cepheid lifetime, meaning the mass changes presented are 
upper limits. 

This estimate suggests that mass loss may be a significant 
contributor but it cannot account for the entire mass discrepancy. 



Table 1. Contribution to the mass discrepancy due to pulsation- 
driven mass loss. 



Mass 


THe 


M 


AM/M 


(M ) 


(Myr) 




(%) 


4 


31.294 


1.4 x 10~ 8 


11 


5 


17.487 


6.4 x 10~ 8 


22 


6 


10.088 


2.6 x 10-* 


4.3 


7 


6.383 


1.3 x 10- 7 


12 


8 


4.867 


3.3 x 10~ 8 


2.0 


9 


3.605 


4.1 x 10- 9 


0.1 



However, detailed calculations are required to explore the feed- 
back of pulsation-driven mass loss on the evolution of Cepheids, 
which in turn, can affect the contribution of mass loss to the mass 
discrepancy. 

3. Stellar evolution models with pulsation-driven 
mass loss 

We computed stellar evolution models for masses M = 
4, 5, 6, 7, 8, and 9 M Q for four different scenarios. The first sce- 
nario is for models with zero convective core overshooting and 
no pulsation-driven mass loss. The second, third, and fourth 
cases include pulsation-driven mass loss and convective core 
overshooting with a c = 0,0.1, and 0.33 5, respect i vely. The 
a c = 0.335 case is based on the results of iBrott et al.l (1201 ll) . In 
the code, convective core overshooting is given by the distance 
that convective cells penetrate above the core on an evolution- 
ary timescale, defined as A = a c Hp, where Hp is the pressure 
scale height and a c a free parameter. Our stellar evolution tracks 
are shown in Fig.Q] The tracks, especially during blue loop evo- 
lution, are very sensitive to the physical processes implemented 
in the mo dels, including both co nvective core overshooting and 
mass loss (Sa lasnich et al.l[l999l) . 

In the models, we assume pulsation-driven mass loss oc- 
curs during the Cepheid stage of evolution only, where the blue 
edge of the Ceph eid instability strip is based on the results of 
Bono et al ] (120001) . and we assume that the red edge is parallel 
to the blue edge with a somewhat arbitrary width. This is be- 
cause the location of the red edge on the HR diagr am is difficult 
to define both observ ationally and theoretically dFernie 1990; 
iFiorentino et al.ll2007l) . However, the chosen location of the red 
edge will not change our main results. Instead of using the 
Leavitt law to determine the pulsat ion period, we now employ 
a period-mass-radius relation from IGieren et al.l (11989k where 
the uncertainty of the predicted period is about 25%. Again, the 
pulsation a mplitudes are computed using the period-amplitude 
relations of iKlagyivik & Szabados! d2009l) . All other fundamen- 
tal parameters describing a Cepheid are taken from the evolution 
models. 

The blue loop evolution of each model differs for each sce- 
nario. When pulsation-driven mass loss is included, the width 
of the blue loops decreases for ma sses M < 1 M Q . This phe - 
nomenon was found previously by iBrunish & Willsonl (|1987), 
who argued that enhanced mass loss can trap a Cepheid in the 
instability strip. 

We can compute the contribution to the mass discrepancy for 
the stellar evolution models. This arises from the contribution to 
the mass discrepancy due to convective core overshoot plus the 
relative change in stellar mass due to mass loss, as shown in 



2 



Neilson et al.: The Cepheid mass discrepancy and pulsation-driven mass loss 



4.5 - 



4 
3.5 

3 
2.5 

o 
_i 

^ 2 
g" 4.5 

4 
3.5 

3 
2.5 





a„ = 





a_ = 0.1 
c i 



a r = 0.335 



4.4 



4.2 



3.8 



3.6 4.4 
LOS T eff 



4.2 



3.8 



3.6 



Fig. 1. Stellar evolution tracks of stars with initial masse s from 4-9 M Q (see labels) from the zero age main sequence until core 
helium exhaustion, computed using the lHeger et al.l ((2000) code for different assumptions of convective core overshooting and mass 
loss during the Cepheid stage of evolution. The dashed lines represent the boundaries of the Cepheid instability strip. 



Table [2] for the three cases with pulsation-driven mass loss in- 
cluded. A Cepheid has a mass discrepancy contri bution o f 2.5% 
for a c = 0.1 and 8.375% for a c = 0.335 from lKelleri d2008). 
For the case of a c - 0.335 with pulsation-driven mass loss, the 
predicted contribution to the m ass discrepan cy, shown in Fig. [2] 
is consistent with the results of iKellerl (12008). This suggests that 
pulsation-driven mass loss in Cepheids is a significant contrib- 
utor towards the mass discrepancy, while allowing for a smaller 
amount of convective core overshooting, _ 0-25 - 0.4 0, as 
opposed to the a c = 0.5 - 1 suggested by iKellej d2008l) . The 
value for a c is also consistent with const raints on convective 
core overshooting in other types of sta rs (Clausen et al. 2010; 
Lovekin & Goupilll2010t ISand berg La cy et al.ll2010t iBrott et al. 
201 lh . 



Most evolution models appear to spend most of the Cepheid 
lifetime near the assumed red edge of the instability strip. At 
this location on the Hertzsprung-Russell diagram, the period- 
amplitude relation overestimates the pulsation amplitude. The 
relation has an intrinsic dispersion caused by the range of pe- 
riods and amplitudes a Cepheid may have as it crosses the in- 
stability strip for a given stellar luminosity. Thus the amplitude 
is an averaged value for a given luminosity. On the other hand, 
nonlinear pulsation models suggest that pulsat ion amplitudes ar e 
largest near the middle of the instability strip dBono et al.l l2000). 



iNeilson & Lester! (|2008) argued that the pulsation-driven mass- 
loss rate increases with decreasing pulsation period or increasing 
pulsation amplitudes, thus suggesting that mass-loss rates at this 
location are also overestimated. Our simplified approach there- 
fore seems to lead to a situation where the amount of total mass 
loss and potentially the width of a Cepheid blue loop depend 
on the assumed location of the red edge of the instability strip. 
However, we argue that this is not the case. To show this, we 
computed a 6 M stellar evolution model with a c = 0.335, as- 
suming that the instability strip is about half as wide as shown in 
Fig-Q] The contribution to the mass discrepancy, in this case, is 
approximately the same as before. This is because a star spends 
most of its Cepheid lifetime at the tip of the blue loop, which is 
close to the helium burning timescale, and the mass-loss rate at 
the tip of the blue loop is approximately the same in both cases. 



ICassisi & Salarisl (1201 lh find that stellar evolution models 
can match the mass of the Large Magellanic Cloud Cepheid 
OGLE-LMC-CEP0227, as measured as part of an eclipsing 
binary system, if one includes convective core overshooting 
with a c — 0.2 u sing the definition of the mass discrepancy 
from lKelled(l2008h . This suggests that the mass discrepancy for 
this particular Cepheid is 5%, contrary to the conclusions of 
ICassisi & Salarisl d201 lh who assume that evolution models with 
moderate convective core overshooting are "standard" models. 
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Table 2. Predicted Cepheid mass discrepancy due to pulsation- 
driven mass loss and different amounts of convective core over- 
shooting as functions of initial mass, M;. 







mass discrepancy (AM/M) 


{M B ) 


a c = 


a c = 0.1 


a c = 0.335 


4 


6.25% 


6.00% (3.50%) 


8.72% (0.34%) 


5 


0.20% 


6.10% (3.60%) 


15.98% (7.60%) 


6 


3.00% 


8.30% (4.80%) 


18.54% (10.16%) 


7 


5.57% 


9.07% (6.57%) 


13.52% (5.14%) 


8 


1.75% 


4.50% (2.00%) 


16.24% (7.86%) 


9 


1.67% 


4.40% (1.90%) 


15.00% (7.62%) 



Notes. Values in parenthesis denote the change in mass due to mass 
loss. For the a c = case the change of mass loss is the predicted mass 
discrepancy. 



(2008) prescription explains a significant portion of the Cepheid 
mass discrepancy but not the entire measured discrepancy. The 
remaining discrepancy can be explained by convective core over- 
shooting, with a value of a c that is consistent with measurements 
in eclipsing binary stars, /3 Cephei stars, and early B-type main 
sequence stars. 

The structure and the width of the model blue loops are also 
found to be affected by the amount of Cepheid mass loss. For 
models with higher values of a c , the mass-loss rates appear to 
be high enough to affect the width of the blue loop, which in 
turn increases the contribution of mass loss to the Cepheid mass 
discrepancy. It is this apparent trapping that causes the differ- 
ences between the predicted mass loss contribution towards the 
mass discrepancy shown in Sects. 2 and 3. 

While our results suggest a possible resolution to the mass 
discrepancy, more observations of dynamic masses of Cepheids 
are needed to constrain this theory. 
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5 6 7 
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Fig. 2. Predicted mass difference due to a c = 0.335 and 
pulsation-driven mass loss as a function of initial stellar mass. 
The horizontal lines and gray regio n represent the average mass 
discrepancy found bv lKellerl (120081) . 



We also note that the result of a c = 0.2 applies for only one low- 
mass Ceph eid and may n ot apply for the higher mass Cepheids 
studied by iKellerl (120081) . This result is similar to the results 
shown in Table [2] for our 4 M models, where the mass dis- 
crepancy is found to vary from about 6% to 8%. While this is 
currently the most precise Cepheid mass known, it is still only 
one Cepheid, and more precise dynamic masses are needed to 
constrain the Cepheid mass discrepancy. 

The results for the contribution towards the mass discrepancy 
found in this section differ from those in Sect. 2. The contribu- 
tion from mass loss tends to increase with increasing values of a c 
because the mass-loss rate depends on the ratio of the luminos- 
ity to mass which increases with a c . Furthermore, an increased 
luminosity means a largerradius, hence a longer predicted pulsa- 
tion period. A longer predicted period suggests larger pulsation 
amplitudes, which again, results in higher mass-loss rates. The 
results in Sect. 2 do not include overshooting, so they miss any 
feedback from overshooting onto the mass loss. 



4. Summary 

We have shown that pulsation-driven mass loss during the 
Cepheid stage of stellar evolution based on the lNeilson & Lester! 
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